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Open and hidden charm production in pA collisions 
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Abstract. We discuss production of charmed mesons and J/\P in p(d)A collisions at high energies. We 
argue that when the saturation scale Q„ characterizing the parton density in a nucleus exceeds the quark 
mass m the naive perturbation theory breaks down. Consequently, we calculate a process of heavy quark 
production in both open and hidden channels in the framework of the parton saturation model (Color 
Glass Condensate). We demonstrate that at RHIC such description is in agreement with experimental 
data on charm production. 

PACS. 24.85.+P Quarks, gluons, and QCD in nuclei and nuclear processes - 14.65.Dw Charmed quarks 



1 Coherent charm production at high energies 

Production of a pair of a quark q and antiqiuark q at 
high energies is characterized by two time scales: produc- 
tion time rp and interaction time r; n t. In pA collisions 
in the center-of-mass frame of the qq pair the production 
time is rp ~ 1/(2to), where to is a quark's mass. In the 
nucleus rest frame, this time is Lorentz time-dilated by 
E g /(2m) where E g is energy of a gluon in a proton from 
which the qq originates. In terms of the Bjorken variable 
x 2 = (mT/\fs)e~ v the production time is rp ~ l/(2Mx2), 
where M is a nucleon mass. At RHIC this corresponds to 
the production time rp ~ 15e v fm. On the other hand, the 
interaction time is r; n t ~ Ra — 7 fm. Therefore, we con- 
clude that in the kinematic region y > at RHIC charm 
is coherently produced by the nuclear color field. 

In high energy QCD coherent color field E of a heavy 
nucleus is described by the Color Glass Condensate 00 
131015,6 ,7,8 ,9 . Being a classical field it scales as ~ 1/g 
with the coupling. Also, the property of geometric scaling 
implies that the only dimensional parameter in the clas- 
sical regime of QCD is the saturation scale Q s . Thus, E ~ 
Ql/g. The saturation scale increases as Q s ~ ^4 1 / 6 e A f/ 2 s A / 4 
with atomic number and energy. When energy accumu- 
lated by the color field at the quark's Compton wavelength 
gEX becomes larger than its rest energy m, the color field 
starts coherent production of qq pairs out of the Dirac sea 
with the probability |l()j : 



(1) 
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Thus, coherent production of heavy quarks by the Color 
Glass Condensate is effective when Q s > m. Analysis 
of a large set of experimental data on the high energy 
hadronic reactions suggest that the saturation scale is 
Q s w lAe 0A5v GeV in the center of the Gold nucleus. 
The charmed quark mass is to c w 1.3 GeV. Therefore, the 
coherent production mechanism switches on at y > and 
becomes dominant at y ~ 2 — 3 in agreement with our 
previous estimate. 

In the next two sections I am going to discuss the phe- 
nomenological implications of the coherent charm produc- 
tion in open and hidden channels. 



2 Open charm production 

In the kinematic region we are interested in rp 3> T- mt the 
quark production can be written in a factorized form as a 
convolution of the valence quark and gluon wave functions 
with the rescattering factor. In a quasi-classical approxi- 
mation the differential cross section for the quark produc- 
tion is pA collisions reads |12II13II14| 
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where I assumed for simplicity that the dominant con- 
tribution comes from interaction of the qq pair with the 
target, while rescatterings of the gluon and the valence 
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Fig. 1. Charmed meson production in dA collisions at yfs = 
200 GeV and y = at RHIC \T7\. Solid line: PYTHIA event 
generator neglecting effect of coherence, dashed line: prediction 
of a model based on gluon saturation/Color Glass Condensate 
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quark are neglected. However in general, they must be 
taken into account as well. In Eq. J5J I used the following 
notations: 

a s C F i a % + (1 - a )%o) ■ i a V + (1 - «)xo) 



7r 2 (ax + (1 — a)x Q ) 2 (ay + (1 — a)x ) 
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Fig. 2. Charmed meson yield: one-particle inclusive cross sec- 
tion integrated over the experimental acceptance region in pr 
and averaged over small bins in rapidity y. 



At high energies gluon density in a nucleus of atomic 
number A is A 1 / 3 times larger than one in a proton. There- 
fore, it makes sense to study the centrality dependence of 
charmed meson yield in pA collisions. The centrality de- 
pendence is usually expressed in terms of the "number of 
binary collisions" iV co ii experienced by the colliding sys- 
tems of nucleons at given impact parameter. The result 
of the calculation is presented in Figure I2|ll| . One ob- 
serves a significant coherence effect in the charmed meson 
production at forward rapidities. Similar result was also 
obtained for the charmed meson production in AA colli- 
sions although in that case the charmed meson yield 
is additionally suppressed by absorption in a hot nuclear 
medium. 



Ki(\x - x \ m) Ki(\y - x \ m) 



Ti" V m-x \\y~x \ 
x[a 2 + (l-a) 2 ]+ K (\x- x \m)K (\y- x \m)^j ,(4) 

where Q s is the saturation scale, m is a quark mass. 

The effect of the high energy evolution of the dipole- 
nucleus scattering amplitude can be taken into account 
by evolving the Glauber-Mueller scattering amplitude 1 — 
e~3- ®° according to the Kovchegov equation |15Ulfi| . In 
practice, since the exact analytical solution to the Kovchegov 
equation is not known, one uses models which are sup- 
posed to satisfy the main properties of that equation. Fol- 
lowing one such model suggested in ^1] we predicted the 
spectrum of charmed mesons in dA collisions, see Figure^ 
Our prediction |1 lj describes the experimental data well 
at transverse momenta px > 1 GeV while the PYTHIA 
event generator based on incoherent charm production 
fails. Moreover, fit of the function A(l +pt/po) ~" to the 
experimental data reveals that the intrinsic transverse mo- 
mentum po 1.32 GeV ^2] is quite close to the saturation 
scale is Q s w 1.4 GeV which indicates consistency of our 
approach. 



3 J/\P production 

In addition to the two time scales characterizing the qq 
pair production rp and Tint, the hidden charm production 
involves another important times scale Tp which is the 
time over which the charmonium bound state is formed. 
In the rest frame of the produced particle this time is of 
the order of inverse binding energy 2 / (M^i — ) (we con- 
centrate on J/SP production). In the nucleus rest frame in 
the RHIC kinematic regime Tp ~ 42 e v fm 3> Tp . We can 
therefore distinguish several kinematic regimes character- 
ized by the relations between the time scales rp, Tj nt and 

Tp. 

At y > the cc pair is produced coherently on a whole 
nucleus and J/H? pair is formed outside it. In that case we 
can calculate the J/'P production cross section as a convo- 
lution of four time-separated amplitudes: gluon emission 
off a valence quark of a proton, splitting of the gluon into 
cc pair, coherent interaction of a cc dipole of a given trans- 
verse size with a nucleus and subsequent formation of the 
J/ty wave function. In the quasi-classical approximation 
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the final result is 1181: 
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where F ee — 5.26 KeV is the leptonic width of J/<P. In 
derivation of (J5J we used the non-relativistic approxima- 
tion to the J/W wave function; the scattering amplitude of 
the cc pair is calculated in large N c approximation; valence 
quarks and the intermediate gluon are treated as specta- 
tors and, finally, there are parametrically small corrections 
due to contributions of the real part of the amplitude and 
the off-diagonal matrix elements. 

To study the nuclear effect in inclusive observables one 
usually defines the nuclear modification factor 



RpA 



daP A /dy 
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In the forward kinematic region Q s 3> M,p. In that case 
10 takes form 



dgP A 
dy 
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which implies the following behavior of the nuclear modi- 
fication factor 



R p a(J/&) 



(8) 



It gets suppressed both as a function of energy/rapidity 
and centrality. At forward rapidities coherence in the nu- 
clear wave function has similar effect on J/W production 
as in the open charm case. 

Unlike in the open charm production case we expect 
that at rapidities y <~ the J/'P inclusive cross section 
in pA collisions is slightly enhanced as compared to that 
in pp collisions (scaled by corresponding ^4). The reason 
is that when rp ~ 2 — 3 fm, J is produced on a nucleus 
by exchange of two gluons with a different nucleons, while 
the same process in pp collisions goes with one gluon ex- 
changed and another one emitted. Additional exchanged 
gluon brings in an additional A 1 / 3, enhancement of the 
scattering amplitude. However, at forward rapidities this 
effect is screened by the multiple gluon exchange which 
leads to suppression of nuclear modification factor as we 
discussed already above. It is important to note, that al- 
though in the kinematic region of interest the production 
time is smaller than the interaction time, still the forma- 
tion time is large Tp > rj nt . Due to color transparency 
the effect of absorption in a cold nuclear medium is only a 
small correction (of the order of a few percents). The effect 
of enhancement occurs in the kinematic region 3> Q s 
where we can expand the exponent in JSJ and derive 



R pA (J/V) = A 1/3 
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Fig. 3. Nuclear modification factor as a function of centrality 
at different rapidities. Solid lines: numerical calculation |18) . 
Data points (preliminary) are from |19| . 



which means that the J/W production is enhanced at slightly 
backward rapidities at RHIC and is stronger in central 
events than for peripheral (since otherwise Q s exceeds 
Afy) . 

It should be emphasized that the prediction of en- 
hancement of the J/<P production at slightly backward ra- 
pidities is a model independent statement. However, since 
in that kinematic region there is no clear separation of 
the time scales our approach (see (0) breaks down. It 
certainly cannot be applied at y < — 1. 

As in the case effect of open charm production the high 
energy evolution can be taken into account by evolving the 
quasi-classical scattering amplitudes with the Kovchegov 
evolution equation. In the Figure we present the result 
of our calculation 18 . 

Another method to express the nuclear effect in the 
total cross sections is to define the variable a such that 
tjpA = A a a pp . It easy to express a it terms of the nuclear 
modification factor R p a, ®. In the Figure 0] we show re- 
sult of our calculation together with the experimental data 
at different energies. Note, that since at SPS and Teva- 
tron energies production time is smaller than the inter- 
action time, we introduced the nuclear absorption factor 
SV = 0.6 EDI. 



4 Conclusions 

We argued that if the quasi-classical coherent color field 
of the Color Glass Condensate is so strong that Q s 3> m, 
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then the production pattern of heavy quarks with mass m 
is the same as the light quarks. At RHIC at the center- 
of-mass energy of y/s = 200 GeV the saturation scale be- 
comes much larger than the charmed quark mass at for- 
ward rapidities y ~ 2—3. Therefore, we expect suppression 
of open charm production in pA and AA as compared to 
pp collisions in much the same as the lighter quarks and 
gluons are suppressed. 
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Although suppression of the nuclear modification fac- 
tor for J/\P production at forward rapidities at RHIC 
follows directly from the corresponding suppression of cc 
production, the detailed mechanism which triggers onset 
of saturation in that case is different due to a specific 
global structure of the J/H? wave function. This mecha- 
nism might manifests itself in a narrow kinematic region 
at which rp ~ 2 — 3 fm (about y <~ 0) making possible 
J/SP production via double gluon exchange. Experimen- 
tally, it will come about as an enhancement of the nuclear 
modification factor in peripheral and semi-central events. 

Coherence effects due to saturation of the nucleus wave 
function will be even stronger at LHC. Since the central 
rapidity interval is shifted by ln(5.5/0.2) = 3.3 with re- 
spect to that of RHIC, we expect that the Color Glass 
Condensate will have a dramatic effect on particle pro- 
duction in general, and charm in particular, at LHC in 
pA and A A and collisions at y > — 3, and even on pp ones 
at somewhat higher rapidities. 



